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ABSTRACT 


A computer program, with the code name TRUMP, was developed in 1965 
by the Lawrence Livermore Laboratory to solve a variety of transient 
and steady~state conduction heat flow problems in simple or complicated 
geometric configurations. The purpose of this investigation was to 
adapt "TRUMP'' to the Naval Postgraduate School IBM/360 Model 67 com- 
puter system. Several heat conduction problems are solved with the 
adapted version of TRUMP and the results compare closely with the 
analytical solutions to these problems. The example problem inputs 
may be used as guides in preparing input data for future problems. 


Some suggestions are also given for further development of the program. 
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L. INTRODUCTION 


A computer program with the code name TRUMP, was developed in 1965 
by Edwards | Reference Hl of the Lawrence Livermore Laboratory (LLL) to 
solve a wide variety of non-linear, conduction heat transfer problems. 
TRUMP has been extensively used and continually improved since its 
initial development and has been successfully adapted to a number of 
computer systems including the CDC/6600 (LLL), IBM/7094 (Bell Labs), 
UNIVAC/1108 (Battelle Northwest, Gulf General Atomics, NBS), and the 
IBM/360 (B.F. Goodrich, NASA Lewis Research Center, Bendix Electro- 
dynamics, Cockerill-Ougee-Providence), to name a few. The objective 
of this thesis project was to adapt TRUMP to the Naval Postgraduate 
School (NPS) IBM/360 Model 67 computer system. 

The B.F. Goodrich Corporation supplied a listing of a version of the 
TRUMP program that they had converted for use in their IBM/360 system. 
Cards were punched from this listing and a version of TRUMP was thus 
secured. 

Mie &.r. Goodrich version was not complete. In addition to errors 
that resulted from punching 3,350 cards, the program also had an align- 
ment error in one subroutine of the program. Also, several options 
that are available in the complete LLL version of TRUMP were missing 
@oeene B.F, Goodrich version. 

With a great deal of assistance from the NPS Computer Facility staff, 
the majority of these problems have been corrected and a running version 
of TRUMP is now available at NPS. 

No, attempt 1S madé in this thesis to give a complete description of 


the TRUMP program. Rather, limitations of the NPS IBM/360 version as 





compared to the LLL CDC/6600 version are discussed. In addition, several 
example problems are solved and may be used as guides in preparing input 
data for future problems. For a complete description of the program, 


reference should be made to the TRUMP users manual [Reference i} . 





If, DESCRIPTION OF THE PROGRAM 


A. GENERAL 

The program solves the general, non-linear, parabolic partial dif- 
ferential heat equation. Transient and steady-state heat conduction 
problems in multidimensions may be solved. [Initial conditions, material 
properties, source and sink strength and boundary conditions may vary 
Moen time, spatial position and/or temperature. 

In order to solve problems with TRUMP, the user specifies input 
fata sin "BLOCK form to 12 input data blocks. A description of each 
of these blocks and their use will be given in Section B 2. 

Following the definition of the problem a model must be constructed. 
Complex geometric shapes are divided into regions and if there is sym- 
metry, the symmetry planes are represented by insulated surfaces in 
order to simplify the model. Regions may have different materials with 
Specified properties. 

The regions may be subdivided into volume elements or nodes having 
arbitrarily any size and shape. Nodes or volume elements must have a 
representative nodal point whose location may be anywhere in the node 
or on the surface of the node. However, in transient problems, for 
maximum accuracy, node shapes and nodal point locations are chosen so 
that the lines joining the nodal points of connected nodes are perpen- 
dicularly bisected by the connected area. In steady-state problems, 
due to the lack of heat capacity effects, the nodal points may be lo- 
eated anywhere in the node without loss of accuracy. A variety of 
boundary conditions may be given to the surface nodes. Constant temper- 


ature, variable temperature, constant heat flux, variable heat flux, 
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forced convection conditions, free convection conditions, radiation 
boundary conditions or combinations of the above may be treated. 

Instead of having very fine zoning of the regions, rather crude 
zoning will give the solution with less effort. Later, the solution 
may be obtained for a finer subdivision of the system. Also, with 
the crude zoning one can easily understand which regions need finer 
zoning. 

The accuracy of the given problem depends on how well the problem 
can be modeled, and the fineness of the spatial subdivisions and time 
steps used in the calculations. 

Allowable problem size limits for the NPS IBM/360 version of TRUMP 


are summarized in Table l. 


edi 





TABLE I 


Problem Size Limits 


Ltem 


Materials 

Reactants 

Nodes 

Internal thermal connections 
External thermal connections 
Boundary nodes 

Heat generation tables 
Initial conditions 

Mass flow connections 

Temp. versus time plot nodes 
Temp. versus plot times 
Remotely dependent properties 


Table lengths 


Data Block 


2 


10 


Lis 


rz 


Pree On 13. Ome 


eZ 


Parameter 


No 
M3 
My 
Ms 
Me 
M7 
Me 
My, 
AIG 
STL 


Ml 


Mio 


Mo 


Size 


ibs: 


BOS. 
950 
60 


20 


50 
50 


10 


nS, 


12 





In Penecralwedilyse@omsistent Set OF units may~be used for the input 
data; output quantities are then given in the same unit system. 

Each material, system node, and boundary node can be identified 
by a non-zero integer number fitting within 5-column fields on the 
data cards. These identifications can be in any sequence. [In order 
to reduce the number of input cards arithmetic sequencing is suggested 
so that the nodes having identical descriptions can be submitted as 
maput on only one card. 

For non-symmetric, symmetric with respect to an axis, or to a 
center, the nodal volumes and the thermal connection areas can be cal- 
euilated by a control value in the input data. 

Initial conditions may be specified in BLOCK 1 of the input data. 
These values will be assigned to all nodes described in BLOCK 4, unless 
other values are specified for individual nodes in BLOCK 9. 

Each internal thermal connection between nodes must be described 
in the input data by specifying the two node identification numbers, 
two connector lengths, and two interface dimensional factors. The 
selection of connector lengths depends on the location of the nodal 
point within each node. These locations may be arbitrary, but in 
transient problems should usually be at the geometric centers of the 
nodes. If the amount of the heat flow depends on surface temperatures, 
then they can be located on the surface. 

A system that exchanges heat with its surrounding must have external 
temperatures specified along its boundary. [In the input data, boundary 
node identification numbers and their temperatures, which may be con- 
stant or a function of time, must be given. The term "boundary node" is 


identical in meaning to an external temperature. 
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Surface temperatures can be specified by using boundary nodes with 
specified temperatures and connecting the surface nodes to a boundary 
node with a very large convection coefficient. In order to reach a 
high degree of accuracy, surface nodes should be selected as either 
zero-volume nodes, or very thin nodes. When a large number of surface 
nodes are connected to the same boundary node and the heat transfer 
coefficients are not tabulated, the number of external connections can 
be reduced by replacing them with connections to a single zero-volume 
node that is connected to the boundary node with a very large convection 
eeerficient. 

In order to generate steady-state solutions, TRUMP must first go 
through a set of transient calculations. The steady-state solution is 
obtained by allowing enough problem time to take the transient to its 
conclusion. However, if only steady-state results are desired, the 
control on the average temperature change during a time step (TVARY) 
can be made very large. This allows the program to neglect transient 
stability limitations on its choice of time steps resulting in a mini- 


mum expenditure of computer time on transient calculations. 


B. USE OF COMPUTER 
1. General 
The data deck must include a problem name card, any number of 
Diock number cards with their input cards, and a data end card. The 
first card of each datadeck must be a problem name card. The symbol "*"! 


feo De 10 column 1, any identification of the problem in columns 2 


eacoughn /l and 73 through 8&0. 
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In the program there are 12 allowed input data blocks, with 
block numbers from 1 through 12. Each block must begin with a block 
mumber card, and except for BLOCKS 1 and 11, must end with a blank 
card. The block number card must begin with the word "BLOCK" in 
columns 1 through 5, and the block number, from Ol through 12, in 
columns 6 and /7. 
Data BLOCK 1 must follow the problem name card; all other data 
blocks may be placed in any order. 
ine Last card of each data deck must be a data end card, with 
the word "ENDED" in columns 1 through 5, and "-1'' in columns 6 and 7. 
Baad. t Lona lecasas COntaining comments without "=" in column 
may precede the name card and will appear on the printout. Any desired 
block description in columns 9 through 80 may be made on the block 
Member card, This will also appear on the printout. 
Z- Block Litem Descriptions 
Basically, one should refer to the TRUMP users manual [Reference 1 
for format details and the original descriptions of the input values. 
In this section the main difference between the IBM/360, Model 67 version 
of TRUMP and the LLL, CDC/6600 version are discussed. 
In BLOCK 1; 
MCYC Maximum allowed number of time steps. If zero, neg- 
ative, or unspecified will be 30,000. 
MSEC Maximum allowed machine time, in seconds. If zero, 
Nesative, Of wnspeciried will be 30,000. 
NPUNCH Indicates that the final values of temperatures of nodes, 
Wemolc MACE LONS O01 Leactants, and, Constance Neat -genel-> 


ation rates of nodes will be punched. Since punched 


es. 











output is not obtained in this version of TRUMP, the 
value of NPUNCH should not be used. 

PREECE, Number of time steps between data output on unit ITAPE 
of problem time and table of node numbers and node 
Lempeotuunec meLieaadd ClO tO OUuLpUut on first, second. 
and last time steps. Only used in TRUMP/6600/DS version. 
Pits eOne TO eel ceMO ts TI Se) 1 iethoS Vers] on: 

ITAPE Output unit on which output data will be written as 
deseribed ine IRi1E.@ Eits Input value is tot an use 
for IBM/360 Model 67 version of TRUMP. 

TIMEP Problem time interval between data output, in addition 
to output on first, second, and last time steps. This 
input control value is also not in use for this version 
Of TRUMP. 

In addition to the above differences, the data transfer option 

is not in use in this versim of TRUMP, 
3. Numerical Output 
The amount of output data is controlled by the value of KDATA, 
However, for every problem, all input data and several other quantities 
such as node volumes and connection areas are written out. The results 
of the first, second, and last time steps are always written out. In 
addition, output data are written out at problem time step intervals 
of IPRINT. For the notations used on the printouts one should refer 
to TRUMP users manual [Reference 1 
4. Evaluation of Results 
Any desired end to the problem may be made with the specified 


input values in BLOCK Ll. These may be maximum problem time, TIMAX, 
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maximum temperature, TMAX, minimum temperature, TMIN, maximum number 
of time steps, MCYC, or maximum number of seconds of machine time, 

MSEC. On the printouts, the value of the Sebi onceut sentinel, KWIT, 
will appear. Table II shows the values of KWIT, and the descriptions 


of these values. 


ly 





KWIT 


10 
Ak 


2 


©AbLE LL 


Values of KWIT 


Cause of Problem Ending 


Tier LOu len same ota t eeached= i iMax 
A temperature exceeding TMAX 

A temperature less than TMIN 
Steady-state 

Data 

Relates to the tape usage 

The number of time steps reached MCYC 
Machine time reached MSEC 

BLOCK 2 or BLOCK 4 missing 
Convergence failure 

Items exceeded the size limits 


Table lengths exceeded the size 
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DESCRIPTIONS OF SUBROUTINES 


THERM 


TALLY 


CHEM 


SPECK 


FINK 


GEN 


SURE 


FLOW 


Sink 


BEEK Z 


SD eOUtLiemLOm materia | propertics and node descrip- 
tions. BLOCKS 2, 4, and 12 are included. 

SUDROW Ene LOL tnittialization, torcaling, and checking 
results, converting regular to special nodes, testing 
for ending the problem. BLOCKS 1 and 9 are included. 
Subroutine for chemical reaction. BLOCK 3 is included. 
Subroutine for heat flow between special nodes and 
other nodes or boundary nodes. JIterates to solve 

set of implicit difference equations when special 
nodes connected to each other. 

Subroutine for heat flow between nodes by conduction, 
convection, and radiation. Treats all nodes as regu- 
lar nodes. BLOCK 5 is included. 

Subroutine for internal heat generation. BLOCK 8 

is included. 

Subroutine for heat flow between surface nodes, bound- 
ary nodes by radiation, free and forced convection 
treats all surface nodes as regular nodes. BLOCKS 6 
and 7 are included. 

Subroutine for mass flow between nodes. Treats all 
nodes as regular nodes. BLOCK 10 is included. 
Subroutine for finding an identification number of 

a material, reactant, node or boundary node in a 

block item list. 

Subroutine for finding the identification numbers 


OL a pair of nodes in a block item List. 


Ie 





oo 
= 
> i Saas aS 
@ 
<p tis = 
= a) 
a 





CLOCK, CLOCKL 


PLOT 


PATCH 


REFER 


SRCON 


SUDEOULIMeSsm@Lor initialization before first time 
step. 

Subroutine for making CRT plots and tables of temp- 
erature versus node location, temperature versus 
time, and node location versus time. BLOCK 11 is 
included. 

Subroutine for converting a number with a decimal 
point to floating number, otherwise substituting 

a specified floating point number. 

Subroutine for finding array subscripts of materials, 
reactants, nodes or boundary nodes referred to by 
identification number in another data block, and 
Writing out diagnostic statements whenever one can 
not be found. 

Subroutine which is related to PATCH. Converts a 
number having decimal point to a floating point 


number. 
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III. EXAMPLE PROBLEMS 


Pe oLAB-MELT PROBLEM 

A semi-infinite slab (-co? x co; ~ coc y Coco, 0 <a Koo), 
has thermal conductivity 1.0 cal/sec-°C, density 1.0 e/em>, heat capacity 
0) cal/sec- C, a melting point of 50.0°C, and a latent heat of fusion 
mme25.0 cal/e. The initial temperature is 100.0°C, and for times iC. 
the surface, z = 0, is maintained at 0.0°C. TRUMP is used to calculate 
the temperature distribution in the slab, and the motion of the solid- 
liquid interface, Zint, for times up to 200 sec. The results are com- 
pared with the exact solution. 


1. Analytic Solution 


From Reference 3. analytic results are: 


Z 
ty eZ VGeciiGa amy FOruO Mh SONG 
OE 
Z . 
Oe oO vert (=a econ CT C 1000C 


24t 
Zint = 0.835-1t fOr =O 


eee eunp solution 

The input data for this problem are shown in Figure 1. In order 
bemmave Output at every 10 cycles, in addition to output on the first, 
second, and last time steps, IPRINT is selected as 10 in BLOCK |. For 
temperature accuracy, TVARY is 1.0°C and the problem is ended at 200 sec 
Dy setting TIMAX equal to 200.0. TONE is 100.0 since the initial tem- 
perature of the slab is 100.0°C. oweuers Since the zero-volume sur- 
face node, node -l, is at One mtor & > 0 its temperatiire is specified 


mm bLOCK 9 by setting the value of TT equal to 0.0. In BLOCK 2, the 


“al 





material name, AMAT, and material identification number are arbitrarily 
selected at SMELT and 50, respectively. Other slab thermal properties 
such as density, heat capacity, thermal conductivity, temperature at 
which the latent heat is released and the latent heat of fusion are also 
specified in BLOCK 2. In order to give node descriptions in BLOCK 4 


miewslab LS divided into nodes as shown below. 


! 






mm 





BeOS 


A total field depth of 100.0 cm is chosen so that at 200 sec no 
appreciable temperature change will take place at the maximum value 
oes iin this figure, the node identification numbers, and the node 
thicknesses are shown. The nodal points of all nodes are at the node 
ereers. In BLOCK 5, internal connection node numbers, lengths of 

heat conduction paths from the nodal points to the connected interfaces, 
Bede the area of the connected interfaces are defined. In BLOCK 6, 

a zero~volume surface node, node -1, is connected to a boundary node, 


node 2001, with a heat transfer coefficient of 1.0 x 108cal/sec-cm*-°C, 


a2 





Oe 





and the external surface area of the surface node is selected as unity. 
In BLOCK 7/7, a constant boundary node temperature is defined as 0.0°C. 


Computer running time for this problem was 65.52 sec. 
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B. FINITE-SLAB PROBLEM (suggested by Professor Paul Pucci) 

A large, 1.0 inch thick steel plate, initially at 80.0°F, is 
suddenly submerged in a fluid at temperature 680.0°F. TRUMP is used 
to calculate the transient temperature distribution in the plate assuming 
aeconstant surface heat transfer coefficient of 1248.0 BTU/hr-f£t2-°F, 
The thermophysical properties of the steel were taken as: thermal con- 
ductivity 26.0 BTU/hr-ft-°F, density 487.0 Lbm/£t?, and specific heat 
0.133 BTU/1bm-°F. 


iene Analytic Solution 





T -T oD SindpL -(\nL)2Fo 
e Gos (Ane 
os (Ay ~ 





2), 
a A,L + Sind,LCos AL 


lel 
Applying the Biot number, Bi =~-—~= 2.0 to the eigen value 


AnL 
problem of Cot(AnL) = ; the two eigen values AL = 1.0769 and 
ba 
ML = 3.6436 can be found. After substituting the appropriate value of 





Ba 


the Fourier number, Fo = 2 = 230.4t in the above series solution and 
considering only the first two terms, the result becomes: 


Be eee Cos(1.0769%) - 0.11836085e°°°*/?¢ cos (3.64368) 


fee lcump Solution 

The input data for this problem are shown in Figure 3. For unit 
consistency, lengths in inches are converted to feet and applied to the 
problem. Since the plate has symmetry with respect to the mid-plane, 
this plane is taken as an insulated plane in the TRUMP solution. So, 
in BLOCK 1 symmetry type indicator, KD, is applied as 1 for the non- 
symmetric case, and the total thickness of the steel plate is considered 
as 0.5 inches. Temperature scale indicator, KT, is selected as 3 


because of the Fahrenheit scale. For temperature accuracy, the desired 
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maximum temperature change in each time step, TVARY, is selected as 
1.0°F. Maximum problem time, TIMAX, is given as 1.0 hr. TONE is set 
equal to 80.0 making the initial temperature of all nodes equal to 80.0°F. 
In BLOCK 2, material name and material number are given as, FSLAB and 
50, respectively. Also, the steel properties such as density, specific 
heat, thermal conductivity are defined. In BLOCK 4 node numbers, mate- 
rials of nodes and their dimensions are given with respect to the equally 
spaced 10 nodes for the 0.5 inches of thickness of the slab. In BLOCK 6, 
the zero-volume surface node, node -1l, is connected to a boundary node, 
node 2001, with a surface heat transfer coefficient of 1248.0 BTU/hr-£t2-°F. 
In BLOCK 7, a boundary node temperature of 680.0°F is given. The com- 
puter running time for this problem was 24.46 sec. 

In order to see the effect of finer zoning on the results, the 
Same problem was also solved using 20 and 40 nodes. However, the results 


were all within 1°F accuracy and not much improvement was obtained with 


the finer grid sizes. Computer running time, on the other hand, was 
m@eeeased from 24.46 sec. to 53./6 sec. and to 1/2.64 sec. for 10, 20, 


and 40 nodes, respectively. 
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Seer ME-INFINITE SLAB PROBLEM 
m Se€mil-intinite slab is initially at a uniform temperature of 100.0°c. 
At time t = 0, the surface temperature is suddenly lowered and main- 
tained at a temperature of 0.0°C. The slab has a thermal conductivity 
of 1.0 cal/sec cm-°C, density of 1.0 e/em, and a heat capacity of 
1.0 cal/sec-°C. The temperattbe distibution im Ghe slab. and. the tera! 
heat flow from the slab are calculated as a function of time for times 
up to 200 sec using TRUMP, and the results are compared with the ana- 


lytical solutions. 


imeeanatytical solution 


T(x,t) = 100 erf ( 


ies 
ee =) ones | t 


2 2eump Solution 


The input data form is shown in Figure 6. in BLOCK 1) problem 
controls such as problem time, TIMAX = 200.00 sec, maximum temperature 
change in each time step, TVARY = OC. ana the initial temperature of 
all nodes, TONE = 100.0°C are given. In BLOCKS 2, 4, 5, 6, 7 and 9, the 
thermal properties, zoning procedure, identification numbers, internal 
and external connection descriptions, boundary node temperatures, and 
initial conditions are defined identically as in the slab-melt problem. 


Computer running time for this problem was 27.05 sec. 


34 





GVIS ALINIANI-IWAS WOL ViIvd INdNE 9 “OLA 


PG a SAV Wee) a ot Clete et a ie ee Cm eto) ie Ua UNS 


0°0 T I~ 
“SNOILIGNOD WiLINI 6 O07 
0° 0 1002 
“S3uUNLVUadH31 3UON AYvVaNNOG 2 YD071¢g 
ey ie O° 7 Cyt 1o00z 1- 
“SNOILdIYISSG GNY SYS@WAN 3GON NOILOSNNOD IvNYaLx3 9 ND074 
Ga OT: O° Oo teum@ee Ye oc 
O°1 O°T O°4 OG 9€ GE 
O°1 oa ag G22, tree ce Te 
Oo O°] O° O°T T€é O€ 
O°T ce O°T Coie Gale ce ic 
Ot aap oni aa iz ve 
O°7 Cron 6*0 SO i Ae T 
O°T O°T eG i i- 
“SNOIAdIYOSSG GNVY SYSSWNN 3GON NOILOSNNOD IWNYILINI Ss 9079 
oe% Ook O°8 OS eel 7 9€ 
Oe Or" O°4 He 5 te 
O°T oct Ste 0G i 6 id 
Ooi O°T or 5 st 
) —_ 
“SNOISNSWIG GNV *S3dAi *SiV1iusavh “SUuauKIIN 30O0N ¥% 90078 
: ual Oot O°T OS qvis 
"S31LY3dOYd WWE3HL GNV *SINANGdWOD TWITW3HI *SUZEWAN ‘SSWVN IWiuslvH 2 0028 
OreO7 ot : 
“STOYLNGS W3IG04d Tt O07 
| °945 002-0 3WIL °10N 
TLOT/S/OL °GIIOS 3LINISNI-INSS NI SO1g LVER LNBSISNVUL » 


55 





me ee an et g 


- enre+e Se a ana aae ea e e aee eee \° “* 
. oe 
' ' 


or 


70 


60 


Ol 
eo 


) 


SC 
cS 
© 


TEMPSRAT LJ R 
@) 


NO 
= 


10 





QO 


LLG 


4 


. TRUMP SOLUTION 
-ANALYTIC SOLUTION 


ee 


16 a1) 





ay 12 
DiS AGG eae 


TEMPERATURE DISTRIBUTION IN THE SEMI-INFINITE SLAB PROBLEM 


36 





a 





TABLE IV 


HEAT FLOW FROM SEMI-INFINITE SLAB 


Time (sec) Qtot CTrump) Qtot (Analytic) 
ir. 20 =119.42 -126.56 
ol ~146.12 -151°86 
205 -178.93 = 63 e86 
3.85 =217 .69 -221.43 
By eye: -264 .69 -267 .69 
8.21 -320.90 =o 23 oo 

2.06 -389.63 -391.79 
iy /0 -472.82 -474.75 
Ze? | -573.36 =379.00 
G5), 02 -694.20 5095.7. 
pes 1S -841.12 -842 .50 
ole. 7 6 -1019.00 -1020.29 

teins, 93 -1234.40 -1235.71 

175.44 -1493.20 -1494.59 

200 .00 -1594.30 =1995.77 
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D. CYLINDER PROBLEM 

Aelita telly lengeecylinder Of tint cradiuemewithea Unit surtace 
temperature, is embedded in a medium ikh zero initial temperature. 
TRUMP is used to calculate the temperature distribution in the medium. 
All thermal properties for the medium have unit values. Results are 
compared with an analytical solution. 

1. Trump Solution 

The input data are shown in Figure 8. Because of the geomet- 

rical shape of cylinder, the symmetry type indicator, KD, in BLOCK l 
is selected as 2. Maximum problem time, TIMAX, is made 50, and the 
medium initial temperature, TONE, is set to 0.0. Also, for temperature 
accuracy, TVARY, is selected as 0.05. In BLOCK 2, unit properties of 
the medium are defined. Nodal point locations are chosen so as to be 
able to compare the TRUMP solution with the given table in Reference 5 


rather easily. The zoning of the medium is shown in the figure below. 


~ 


, 2 1.6 ae 32 
5 ie te i ae ( =a + si aa pa 
| | 
| | 
34 
12ly Sin | 24 iy) 2 fy] 34 


__— 
\ 
‘ — 
\ 
\ 
\ 


By, 





In BLOCK 4, the arithmetic mean average radii, DRAD, are used to obtain 

correct nodal volumes for cylindrical shaped nodes. Actual node inter- 

faces are used as interface radii, DRAD, in BLOCK 5, to obtain the areas 
of the connected interfaces. In BLOCK 6, zero-volume surface node, 

node -1, is connected to the boundary node with a heat transfer coeffi- 


erent of 1.0 x 108 | and geometric factors for external surface node 


area are defined. Computer running time for this problem was 15.46 sec. 
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TABLE V 


TEMPERATURE DISTRIBUTION FOR CYLINDER (TIME = 50) 


tT (Trump) T (Analytic) 

16a e e.O) 

OF IcL OEI6S 
Oreo O) otsene. 
0.842 0.843 
0.794 0.794 
OVE Yes) SNe 
0.574 O39 74 
Clio ORs i 
Seales : 0.260 
On ey) OL 77 
0.146 0.146 
0.016 OF OLG 
O) ste KOVAL 0.001 
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Pee Sen ERE PROBLEM 

A sphere with a unit radius and a unit surface temperature is 
embedded in a medium with a zero initial temperature. The temperature 
distribution in the medium is calculated using TRUMP assuming all the 
thermal properties have unit values. Results are compared with the ana- 


lytical solution. 


Pe Analytic Solution 


rol 
T(r,t) = erfe (5 


2. Trump Solution 
The input data from the cylinder problem can be used to solve 
this problem if the symmetry type indicator, KD, in BLOCK 1] is changed 
from 2 to 3 to correspond to the centrisymmetric geometry. All other 
input conditions are identical with the cylinder problem. The computer 


time used for this problem was 11.6] sec. 
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TABLE VI 


TEMPERATURE DISTRIBUTION FOR SPHERE (TIME = 50) 


T_ (Trump ) T (Analytic) 
E30 ao 
Grog 0.90 
Ory > O75 
0.64 0.64 
OASYS 0.56 
0.49 0.49 
O2Zo 0.28 
0.14 0.14 
0.08 0.08 
CUS 0.05 
0.04 0.04 
0.00 OZ 00 
0.00 0.00 
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F. TWO DIMENSIONAL PLATE PROBLEM 

A rectangular plate, 40 x 20 cm, has a thermal conductivity of 
400 cal/cm-sec- C, a density of 1.0 eieriee and a heat capacity of 
iO cal/e--C ae LS dbMietetlive ele 0) USC! «Ube cules sor: idee jollaue cies 
maintained at 0.0°C. The upper side of the plate has a temperature 
fistribution that varies sinusoidally with x, starting at O°C, reaching 
a maximum of 100°C at the midpoint and decreasing to O°C at the outer 
edge (see figure below). The steady-state temperature distribution is 
calculated in the plate using TRUMP. Results are compared with the 


analytical solution. 
+ 
| : 


T = 100Sin 4) 


~—+ 

U1 i 

O, EOS aeRO ee 
C) 





ie ee 
Btad fe 40cm 


Meenalytie Solution 
Tt y: 
Sr win ae) 
T(x,y) = 100 eee ee iC) 
Sinh(-~) 


where 0 <x Ke 40 and 0 <y € 20 
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Pike Trump Solution 


The input data are shown in Figure 9. Because the boundary 
conditions of the plate and the plate itself are symmetric with respect 
to the plane at x = 20 cm, this symmetry plane is assumed to be an in- 
sulated surface. The zoning of the plate is shown in Figure 10. In 
BLOCK Porlthe Anput data,  EPRMNT, is selected as99999, so that print- 
outs are obtained only for the first, second, and last time steps. To 
obtain a rapid conclusion to the transient problem, TVARY, is selected 
as 100.0°C. In BLOCK 4, instead of using a very large number of surface 
nodes for sides x = 40 cm and y = 0, surface node -1 and surface node -4 
are used, respectively, in order to reduce the number of the input cards, 
and in BLOCK 6 they are connected to the boundary nodes 2001 and 2004. 
In BLOCK 7, the boundary node temperatures, and in BLOCK 9, the initial 
temperatures of the surface nodes are given. The computer running time 


for this problem was 28.71 sec. 
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IV. CONCLUSIONS 


The adaptation of the computer program "TRUMP" to the Naval Post- 
graduate School IBM/360 Model 67 has been achieved, with success. 
Analytical temperature distributions of various transient and steady- 
state heat conduction problems compared very well with the results 
obtained from TRUMP. The complete adaptation has not yet been estab- 
lished, however. At present, punched output and plotted graphs have 
not been obtained. The computer solutions of the example problems did 
not require the use of subroutines CHEM, GEN, and FLOW. Thus, the 
adaptation of these subroutines is yet to be evaluated. Since, the 
accuracy of the solutions are dependent on the zoning, this effort can 
be reduced if an existing program "FED: A computer program to generate 
geometric input for the heat-transfer code TRUMP" by Dale A. Schauer 
Reference 4 is adapted to the IBM/360 Model 67 version of TRUMP. 


The example problems presented herein will provide a useful reference 


Peoeepreparation of input data for future problems. 
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